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Abstract 

Relativistic energy density functionals (REDF) provide a complete and accurate, global descrip- 
tion of nuclear structure phenomena. A modern semi-empirical functional, adjusted to the nuclear 
matter equation of state and to empirical masses of deformed nuclei, is applied to studies of shapes 
of superheavy nuclei. The theoretical framework is tested in a comparison of calculated masses, 
quadrupole deformations, and potential energy barriers to available data on actinide isotopes. 
Self-consistent mean-field calculations predict a variety of spherical, axial and triaxial shapes of 
long-lived superheavy nuclei, and their alpha-decay energies and half-lives are compared to data. 
A microscopic, REDF-based, quadrupole collective Hamiltonian model is used to study the effect 
of explicit treatment of collective correlations in the calculation of Q a values and half-lives. 

PACS numbers: 21.60.Jz, 21.10.Dr, 21.10.Re, 24.75. +i, 27.90. +b 
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I. INTRODUCTION 



For many years theoretical studies of superheavy nuclei (SHN) were mostly based on the 
traditional macroscopic-microscopic approach [l-5], but since the late 1990s the framework 
of self-consistent mean-field models, based on realistic effective inter-nucleon interactions 
or energy density functionals, has systematically been applied to the structure of SHN 



231 ] . Binding energies, deformations, a-decay energies and half-lives, fission barriers and 
spontaneous-fission half-lives, fission isomers, and single-nucleon shell structure of SHN have 
successfully been described using self-consistent mean-field (SCMF) models based on the 
Gogny effective interaction, the Skyrme energy functional, and relativistic meson-exchange 
effective Lagrangians. 

The advantages of using SCMF models include the intuitive interpretation of results in 
terms of single-particle states and intrinsic shapes, calculations performed in the full model 
space of occupied states, and the universality that enables their applications to all nuclei 
throughout the periodic chart. The latter feature is especially important for extrapolations 
to regions of exotic short-lived nuclei far from stability for which few, if any, data are 
available. In addition, the SCMF approach can be extended beyond the static mean-field 
level to explicitly include collective correlations and, thus, perform detailed calculations of 
excitation spectra and transition rates. 

During the last decade, important experimental results on the mass limit of the nuclear 
chart have been obtained using compound nucleus reactions between the 48 Ca beam and 
actinide targets. A number of isotopes of new elements with the atomic number Z = 
113 — 118 have been dicovered, and new isotopes of Z = 110 and 112 24j-l3l|. The decay 



energies and the resulting half-lives provide evidence of a significant increase of stability 
with increasing neutron number in this region of SHN. Theoretical studies predict that SHN 
in this region should display rapid shape transitions, from prolate, through spherical, to 



oblate-deformed ground states 
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321 ] . These nuclei, therefore, present an 



ideal testing ground for structure models attempting to predict the location of an "island of 
stability" for SHN around N = 184. 

In this work we apply the framework of relativistic energy density functionals (REDF) to 
an illustrative study of shape transitions and shape coexistence in SHN with Z = 110 — 120. 
There are several advantages in using functionals with manifest covariance and, in the con- 



text of this study, the most important is the natural inclusion of the nucleon spin degree 
of freedom, and the resulting nuclear spin-orbit potential which emerges automatically with 
the empirical strength. Our aim is to test the recently introduced functional DD-PC1 33| in 
self-consistent relativistic Hartree-Bogoliubov (RHB) calculations of energy surfaces (axial, 
triaxial, octupole), a-decay energies and half-lives of SHN, in comparison to available data 
and previous theoretical studies. Section [TT] introduces the general framework of REDFs and, 
in particular, the functional DD-PC1 that will be used in illustrative calculations through- 
out this work. The functional is tested in calculations of binding energies, ground-state 
quadrupole deformations, fission barriers, fission isomers, and Q a values for even-even ac- 
tinide nuclei. In Sec. IHII we apply the RHB framework based on the functional DD-PC1 and 
a separable pairing interaction in a description of triaxially deformed shapes and shape tran- 
sitions of even-even superheavy nuclei. The microscopic, REDF-based, quadrupole collective 
Hamiltonian model is used to study the effect of explicit treatment of collective correlations. 
Q a values and half-lives for two chains of odd-even and odd-odd superheavy systems are 
computed using a simple blocking approximation in axially symmetric self-consistent RHB 
calculations. Section [TV] summarizes the results and presents an outlook for future studies. 



II. THE RELATIVISTIC ENERGY DENSITY FUNCTIONAL DD-PC1 

Relativistic energy density functionals (REDF) provide an accurate, reliable, and consis- 
tent description of nuclear structure phenomena. Semi-empirical functionals, adjusted to a 
microscopic nuclear matter equation of state and to bulk properties of finite nuclei, are ap- 
plied to studies of arbitrarily heavy nuclei, exotic nuclei far from stability, and even systems 
at the nucleon drip-lines. REDF-based structure models are being developed that go beyond 
the mean-field approximation, and include collective correlations related to restoration of 
broken symmetries and to fluctuations of collective variables. 

Although it originates in the effective interaction between nucleons, a generic density 
functional is not necessarily related to any given nucleon- nucleon potential and, in fact, 
some of the most successful modern functionals are entirely empirical. Until recently the 
standard procedure of fine-tuning global nuclear density functionals was to perform a least- 
squares adjustment of a small set of free parameters simultaneously to empirical properties 
of symmetric and asymmetric nuclear matter, and to selected ground-state data of about 
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ten spherical closed-shell nuclei. A new generation of semi-microscopic and fully microscopic 
functionals is currently being developed that will, on the one hand, establish a link with the 
underlying theory of strong interactions and, on the other hand, provide accurate predictions 
for a wealth of new data on short-lived nuclei far from stability. To obtain unique param- 
eterizations, these functionals will have to be adjusted to a larger data set of ground-state 



properties, including both spherical and deformed nuclei [34J, |35 |. 

For a relativistic nuclear energy density functional the basic building blocks are densities 
and currents bilinear in the Dirac spinor field ip of the nucleon: ipOrTip , O t G {l,Tj}, 
T G {1, 7 M , 75, 757^, cr^j,}. are the isospin Pauli matrices and Y generically denotes the 
Dirac matrices. The isoscalar and isovector four-currents and scalar density are defined as 
expectation values of the corresponding operators in the nuclear ground state. The nu- 
clear ground-state is determined by the self-consistent solution of relativistic Kohn-Sham 
single-nucleon equations. To derive those equations it is useful to construct an interaction 
Lagrangian with four-fermion (contact) interaction terms in the various isospace-space chan- 
nels: isoscalar-scalar (ipip) 2 , isoscalar-vector (^7 /t ?/>)('?/>7' 1 -?/>), isovector-scalar (iprip) ■ (ipTip), 
isovector- vector [ipf^^ip) ■ (i[)T^ip). A general Lagrangian can be written as a power series in 
the currents ^0 T Tip and their derivatives, with higher-order terms representing in-medium 
many-body correlations. 

In Ref. 33| a Lagrangian was considered that includes second-order interaction terms, 
with many-body correlations encoded in density-dependent strength functions. A set of 10 
constants, that control the strength and density dependence of the interaction Lagrangian, 
was fine-tuned in a multistep parameter fit exclusively to the experimental masses of 64 
axially deformed nuclei in the regions A w 150 — 180 and A w 230 — 250. The resulting 
functional DD-PC1 has been further tested in calculations of binding energies, charge radii, 
deformation parameters, neutron skin thickness, and excitation energies of giant monopole 
and dipole resonances. The corresponding nuclear matter equation of state is characterized 
by the following properties at the saturation point: nucleon density p sa t = 0.152 fm~ 3 , 
volume energy a v = —16.06 MeV, surface energy a s = 17.498 MeV, symmetry energy 
04 = 33 MeV, and the nuclear matter compression modulus K nm = 230 MeV. 

For a quantitative description of open-shell nuclei it is necessary to consider also pairing 
correlations. The relativistic Hartee-Bogoliubov (RHB) framework [36j provides a unified 
description of particle-hole (ph) and particle-particle (pp) correlations on a mean-field level 



by combining two average potentials: the self-consistent mean field that encloses all the 
long range ph correlations, and a pairing field A which sums up the ^-correlations. In this 
work we perform axially-symmetric and triaxial calculations based on the RHB framework 
with the ph effective interaction derived from the DD-PC1 functional. A pairing force 
separable in momentum space: (k\ V ls ° \k') = —Gp(k)p(k') is used in the pp channel. By 
assuming a simple Gaussian ansatz p{k) = e~ a2fc2 , the two parameters G and a were adjusted 
to reproduce the density dependence of the gap at the Fermi surface in nuclear matter, 



calculated with a Gogny force. For the D1S parameterization 37| of the Gogny force the 



38|. When 



following values were determined: G = —728 MeVfm 3 and a = 0.644 fm 
transformed from momentum to coordinate space, the force takes the form: 

V(r lf r 2 , ri, r' 2 ) = G5(R- X) P(r)P(r')± (1 - P ff ) , (1) 

where R = \ [r\ + r%) and r = r\ — r 2 denote the center-of-mass and the relative coordi- 
nates, and P(r) is the Fourier transform of p{k): P{r) = 1/ (Atto 2 ) 3 ^ 2 e^ 7 * 2 / 4 " 2 . The pairing 
force has finite range and, because of the presence of the factor 5(R — R), it preserves 
translational invariance. Even though 5 (R — R') implies that this force is not completely 
separable in coordinate space, the corresponding anti-symmetrized pp matrix elements can 
be represented as a sum of a finite number of separable terms in the basis of a 3D harmonic 
]. The force Eq. ([!]) reproduces pairing properties of spherical and deformed 



oscillator 



nuclei calculated with the original Gogny force, but with the important advantage that the 
computational cost is greatly reduced. 



The Dirac-Hartree-Bogoliubov equations 36J are solved by expanding the nucleon spinors 



in the basis of a 3D harmonic oscillator in Cartesian coordinates. In this way both axial 
and triaxial nuclear shapes can be described. The map of the energy surface as a function 
of the quadrupole deformation is obtained by imposing constraints on the axial and triaxial 
quadrupole moments. The method of quadratic constraint uses an unrestricted variation of 
the function 

(H)+J2 C2,((Q 2 ,)-^) 2 , (2) 

M=0,2 

where (H) is the total energy, and (C^) denotes the expectation value of the mass 
quadrupole operators: 

Q 20 = 2z 2 - x 2 - y 2 and Q 22 = x 2 - y 2 . (3) 
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FIG. 1. (Color online) Absolute deviations of the self-consistent RHB ground-state binding energies 
from the experimental values [40], for the isotopic chains of Th, U, Pu, Cm, Cf, Fm, and No (upper 
panel). In the lower panel the calculated ground-state axial quadrupole moments are shown in 



comparison to data 



41] (open symbols). 



q2 fl is the constrained value of the multipole moment, and the corresponding stiffness 
constant. 

To illustrate the accuracy of the DD-PC1 functional in the calculation of ground-state 
properties of heavy nuclei, in Fig. [T]we plot the results of self-consistent 3D RHB calculations 
for several isotopic chains in the actinide region: Th, U, Pu, Cm, Cf, Fm, and No. The 

n 

deviations of the calculated binding energies from data [4CJ show an excellent agreement 
between theory and experiment: the absolute difference between calculated and experimental 
binding energies is less than 1 MeV in all cases. An important results is also that the 
mass residuals do not display any notable dependence on the mass (neutron) number. The 



calculated ground-state quadrupole Q20 moments are compared with available data |4l| in 
the lower panel of Fig. [TJ One notices that the values predicted by the DD-PC1 functional 
reproduce in detail the isotopic trend of the empirical moments in the Th, U, Pu and Cm 
sequences, and are in very good agreement with the quadrupole moments of the Cf isotopes. 

The "double-humped" fission barriers of actinide nuclei provide an important test for 
nuclear energy density functionals. In the review of self-consistent mean-field models for 



nuclear structure 15[, which also contains an extensive list of references to previous studies 
of fission barriers using mean-field-based models, Bender et al. compared paths in the 
deformation energy landscape of 240 Pu obtained with various Skyrme, Gogny and relativistic 
mean-field (RMF) interactions. In general, relaxing constraints on symmetries lowers the 
fission barriers. The predicted shapes are triaxial and reflection-symmetric at the first 
barrier, and predominantly axial and reflection-asymmetric at the second barrier. The 
systematics of axially symmetric fission barriers in Th, U, Pu, Cm and Cf nuclei, as well 
as for superheavy elements Z = 108 — 120, using several Skyrme and RMF mean-field 



interactions, was investigated in Ref. 42]. The fission barriers of twenty-six even-Z nuclei 



with Z = 90 — 102, up to and beyond the second saddle point, were calculated in Ref. [43 1 
with the constrained Hartree-Fock approach based on the Skyrme effective interaction SkM*. 
The fission barriers of 240 Pu beyond the second saddle point were also explored using the 



axially quadrupole constrained RMF model with the PK1 effective interaction 44j. In a 
very recent optimization of the new Skyrme density functional UNEDF1 {35}], excitation 
energies of fission isomers in 236 > 238 u, 240 Pu, and 242 Cm, were added to the data set used 
;o adjust the parameters of the functional. Compared to the original functional UNEDFO 



34j, the inclusion of the new data allowed an improved description of fission properties of 
actinide nuclei. The effect of triaxial deformation on fission barriers in the actinide region 
was recently also explored in a systematic calculation of Ref. [45], based on the RMF+BCS 
framework. The potential energy surfaces of actinide nuclei in the P22, ^30) deformation 
space (triaxial + octupole) were analyzed in the newest self-consistent mean-field plus BCS 



calculation based on relativistic energy density functionals [46|. This study has shown the 
importance of the simultaneous treatment of triaxial and octupole shapes along the entire 
fission path. 

The fission barriers calculated in the present work are shown in Fig. |2j where we plot the 
potential energy curves of 236 > 238 U, 240 Pu, and 242 Cm, as functions of the axial quadrupole de- 



formation parameter fao ■ The deformation parameters are related to the multipole moments 
by the relation: 

An 

= zXj&iQ*) • ( 4 ) 

To be able to analyze the outer barrier heights considering reflection-asymmetric (octupole) 
shapes, the results displayed in this figure have been obtained in a self-consistent RMF plus 
BCS calculation that includes either triaxial shapes, or axially symmetric but reflection- 
asymmetric shapes. The interaction in the particle-hole channel is determined by the rela- 
tivistic functional DD-PC1, and a density-independent 5-force is the the effective interaction 



• 3, 



in the particle-particle channel. The pairing strength constants V n and V p are from Ref. 
where they were adjusted, together with the parameters of the relativistic functional PC-F1, 
to ground-state observables (binding energies, charge and diffraction radii, surface thickness 
and pairing gaps) of spherical nuclei, with pairing correlations treated in the BCS approx- 
imation. In many cases the functionals PC-F1 (471 and DD-PC1 predict similar results for 
ground state properties (cf. for instance, Ref. 481]). and reproduce the empirical pairing 
gaps. Thus we assume that, without any further adjustment, the same pairing strength 
parameters can be used in RMF+BCS calculations with the functional DD-PC1. 

The solid (black) curves correspond to binding energies calculated with the constraint on 
the axial quadrupole moment, assuming axial and reflection symmetry. The absolute minima 
of these curves determine the energy scale (zero energy). The dot-dashed (blue) curves 
denote paths of minimal energy in calculations that break axial symmetry with constraints 
on quadrupole axial Q20 an d triaxial Q22 moments. Finally, the dashed (green) curves are 
paths of minimal energy obtained in axially symmetric calculations that break reflection 
symmetry (constraints on the quadrupole moment Q20 and the octupole moment Q30). The 
red squares, lines, and circles denote the experimental values for the inner barrier height, 
the excitation energy of the fission isomer, and the height of the outer barrier, respectively. 
The data are from Ref. {49 1. 

The excitation energies of fission isomers are fairly well reproduced by the axially sym- 
metric and reflection symmetric calculation, but the paths constrained by these symmetries 
overestimate the height of the inner and outer barriers. The inclusion of triaxial shapes 
lowers the inner barrier by ~ 2 MeV, that is, the axially symmetric barriers in the region 
020 ~ 0.5 are bypassed through the triaxial region, bringing the height of the barriers much 
closer to the empirical values. As shown in the figure, the inclusion of octupole shapes 
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FIG. 2. (Color online) Constrained energy curves of 236 > 238 u, 240 Pu, and 242 Cm, as functions 
of the axial quadrupole deformation parameter. Results of self-consistent axially and reflection- 
symmetric, triaxial, and axially reflection-asymmetric RMF+BCS calculations are denoted by solid 
(black), dot-dashed (blue), and dashed (green) curves, respectively. The red squares, lines, and 
circles denote the experimental values for the inner barrier height, the excitation energy of the 



fission isomer, and the height of the outer barrier, respectively. The data are from Ref. 



49] 



(axial, reflection-asymmetric calculations) is essential to reproduce the height of the outer 
barrier in actinide nuclei. A very good agreement with data is obtained by following paths 
through shapes with non-vanishing octupole moments. With the present implementation of 
the model we could not simultaneously calculate both octupole and triaxial shapes, but such 
a study was recently performed in the RMF+BCS framework by Bing-Nan Lu and collabo- 
rators 46j. It was shown that not only the inner barrier, but also the reflection-asymmetric 



outer barrier is lowered by the inclusion of triaxial deformations. This effect is of the order 



> 

CD 




148 150 
N 



FIG. 3. (Color online) Q a values for even-even actinide chains obtained in a self-consistent axially 
symmetric RHB calculations using the functional DD-PC1 and the separable pairing interaction 
Eq. (pQ). The theoretical values (filled symbols) are connected by lines and compared to data (open 



symbols) 



4C] 



of 0.5 - 1 MeV, and it accounts for 10% to 20% of the barrier height. Considering that the 
functional DD-PC1 was adjusted only to the binding energies of the absolute axial minima 
(masses) of deformed nuclei, the results shown in Fig. [2] reproduce the experimental values 
surprisingly well and, therefore, appear to be very promising for its extrapolation to the 
region of superheavy nuclei. We note here that, except for the results of Fig. [2} all the other 
calculations reported in this work have been performed in the RHB framework using the 
functional DD-PC1 and the separable pairing interaction Eq. (pQ). 

Figure [3] illustrates the accuracy of the functional DD-PC1 in the axially symmetric RHB 
calculation of Q a values, that is, energies of a particles emitted by even-even actinide nuclei. 



The calculated values are plotted in comparison to data |40j]. Even in this simple calculation 
that assumes axial symmetry, the model reproduces the empirical trend of Q a values. The 
few cases for which we find a somewhat larger deviation from data most probably point 
to a more complex potential energy surface, possibly including shape coexistence. It is 
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interesting to note that on the quantitative level the theoretical results are very similar 
to those obtained in the self-consistent non-relativistic Hartree-Fock-Bogoliubov calculation 



based on the Skyrme functional SLy4 [13l. Il7|. 



Summarizing this section, it has been shown that self-consistent mean-field calculations 
based on the relativistic energy density functional DD-PC1 predict binding energies, ground- 
state quadrupole deformations, fission barriers, fission isomers, and Q a values for even-even 
actinide nuclei in very good agreement with data. In the next section we apply the RHB 
framework based on the functional DD-PC1 and the separable pairing interaction Eq. ([1]) 
to an illustrative study of deformed shapes and shape transitions of superheavy nuclei. 



III. SHAPE TRANSITIONS IN SUPERHEAVY NUCLEI 



In a very recent study of fission barriers and fission paths in even-even superheavy nuclei 



with Z = 112 — 120 [23j, DD-PC1 was used, together with two other relativistic energy 
density functionals, in a systematic RMF+BCS calculation of potential energy surfaces, 
including triaxial and octupole shapes. It was shown that low-Z and low-N nuclei in this 
region are characterized by axially symmetric inner fission barriers. With the increase of Z 
and/or N, in some of these nuclei several competing fission paths appear in the region of the 
inner barrier. Allowing for triaxial shapes lowers the outer fission barrier by 1.5 — 3 MeV, 
and in many nuclei the lowering induced by triaxiality is even more important than the one 
due to octupole deformation. 

The variation of ground-state shapes is governed by the evolution of the shell struc- 
ture of single-nucleon orbitals. In very heavy deformed nuclei, in particular, the density of 
single-nucleon states close to the Fermi level is rather large, and even small variations in 
the shell structure predicted by different effective interactions can lead to markedly distinct 
equilibrium deformations. To illustrate the rapid change of equilibrium shapes for the heav- 
iest nuclear systems, Figure S] displays the results of self-consistent axially symmetric RHB 
calculations of isotopes in the a-decay chains of two superheavy nuclei: 298 1 20 and 300 1 20. 
The quadrupole energy curves are plotted as functions of the deformation parameter f3 2 o ■ 
Lighter systems around Z = 110 are characterized by well developed prolate minima around 
/?2o ~ 0.2, whereas intermediate nuclei display both oblate and prolate minima at small 
deformation, and the heaviest isotopes appear to be slightly oblate. Another characteristic 
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of the energy curves is the shift of the saddle point to smaller deformations with the increase 
in the mass number, while the barriers become wider. Since the prolate and oblate min- 
ima can be connected through triaxial shapes without a barrier, these energy curves show 
the importance of performing more realistic calculations, including triaxial, and for large 
deformations, octupole shapes. 



168 170 172 174 176 178 180 




FIG. 4. (Color online) Self-consistent RHB axially symmetric energy curves of isotopes in the 
o-decay chains of 298 120 and 300 120, as functions of the quadrupole deformation parameter. 



This is shown in Figs. [5] and EJ where we plot the corresponding triaxial RHB energy 
surfaces in the /3 — 7 plane (0 < 7 < 60°) for isotopes in the a-decay chains of 298 120 and 
300 120, respectively. In both chains the heaviest systems display soft oblate axial shapes 
with minima that extend from the spherical configuration to |/?2o| ~ 0.4 [Z = 120) and 
l&o I ~ 0.3 (Z = 118). We do not have to consider the deep prolate minima at /3 2 o > 0.5 
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FIG. 5. (Color online) Self-consistent RHB triaxial energy maps of the even-even isotopes in the 
a-decay chain of 298 120 in the (3 — 7 plane (0 < 7 < 60°). Energies are normalized with respect to 
the binding energy of the absolute minimum. 
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because the inclusion of reflection asymmetric shape degrees of freedom (octupole deforma- 
tion) drastically reduces or removes completely the outer barrier. A low outer barrier implies 
a high probability for spontaneous fission, such that the prolate superdeformed states are 



not stable against fission 50J. In contrast to the actinides shown in Fig. [2J superheavy nuclei 
are actually characterized by a "single-humped" fission barrier. As already noted, in the 
present implementation of the model triaxial and octupole deformations cannot be taken into 
account simultaneously. Reflection asymmetric shape degrees of freedom, however, play no 
role at small and moderate deformations that characterize ground-state configurations of the 
superheavy systems considered here. The intermediate nuclei with Z = 116 are essentially 
spherical but soft both in /3 and 7, whereas prolate deformed mean-field minima develop 
in the lighter systems with Z — 114, Z — 112 and Z = 110. The predicted evolution of 
shapes is consistent with results obtained u sing the self-consistent Hartree-Fock-Bogoliubov 
framework based on Skyrme functionals {(j. lid u\. 

The two main decay modes in this region are a-emission and spontaneous fission. The 
theoretical a-decay energies, denoted by (blue) diamonds in Fig. [7J are calculated as the 
difference between the absolute minima of the energy maps of the parent and daughter 
nuclei, shown in Figs. |5] and [6j For the mean-field ground state we take the minimum 
with the highest barrier with respect to fission, that is, we do not consider superdeformed 
minima with very low fission barriers. The theoretical Q a values are shown in comparison 



to available data for the a decay properties of superheavy nuclei |24j . The trend of the data 
is obviously reproduced by the calculation, and the largest difference between theoretical 
and experimental values is less than 1 MeV. This is a rather good result, considering that 
equilibrium nuclear shapes change rapidly in the two a-decay chains and that the calculation 
is performed on the mean-field level. In general the level of agreement with experiment is 
similar to the one found in the case of actinide nuclei (cf. Fig. |3]), but not quite the same 
as that obtained using the macro-micro model specially adapted to heaviest nuclei (HN) 



ABES 



-53|. It is interesting, though, that our prediction both for the a-decay energy and 
half-life of the A = 298 isotope of the new element Z = 120 are consistent with those of the 
HN macro- micro model 53 1. 



Alpha-decay half-lives are calculated using a simple five-parameter phenomenological 
Viola-Seaborg-type formula 0, Q. The parameters of this formula were adjusted to experi- 
mental half-lives and Q a values of more than 200 nuclei with Z = 84 — 111 and iV = 128 — 161 
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FIG. 6. (Color online) Same as described in the caption to Fig. [5] but for the the a-decay chain of 
300 120. 



54J. Using the theoretical Q n values plotted in Fig. [7] as input, the resulting half-lives are 



compared to available data 24j in Fig. El A rather good agreement with experiment is 
obtained for both decay chains. 

The theoretical values denoted by (blue) diamonds in Figs. [7J and correspond to tran- 
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FIG. 7. (Color online) Q a values for the a-decay chains of 298 120 (left) and 300 120 (right). The 
theoretical values are calculated as the difference between the mean-field minima of the parent and 
daughter nuclei (blue diamonds), and as the difference between the energies of the + ground states 



of the quadrupole collective Hamiltonian (red circles). The data (squares) are from Ref. 



24]. 
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FIG. 8. (Color online) Half-lives for the a-decay chains of 298 120 (left) and 300 120 (right). The 



54|, using 



theoretical values are calculated from a phenomenological Viola-Seaborg-type formula 5, 
the Q a values from Fig. [71 Diamonds correspond to values of Q a calculated from mean-field RHB 
solutions, whereas circles denote half-lives computed using Q a values determined by the + ground 



states of the quadrupole collective Hamiltonian. The data (squares) are from Ref. 



24]. 



sitions between the self-consistent mean-field minima on the triaxial RHB energy surfaces 
shown in Figs. [5] and EJ Such a calculation does not explicitly take into account collective 
correlations related to symmetry restoration and to fluctuations in the collective coordinates 
and 7. Physical transitions occur, of course, not between mean-field minima but between 
states with definite angular momentum. In cases in which both the initial and final states 
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FIG. 9. (Color online) Q a values (left), and half-lives (right) for the a-decay chain of 287 115, and 
for the nucleus 293 117. 



have similar deformation this will not make a large difference for the Q a values, because 
collective correlations are implicitly taken into account in energy density functionals through 
the adjustment of parameters to ground-state properties (masses). The difference, however, 
can be larger in cases when the equilibrium shapes of the parent and daughter nucleus cor- 
respond to rather different deformations, because the collective correlation energy generally 
increases with deformation. For this reason we have also used a recent implementation of the 
collective Hamiltonian based on relativistic energy density functionals 55], to calculate a- 
transition energies between ground states of even-even nuclei (0 + — > + transitions). Starting 
from self-consistent single-nucleon orbitals, the corresponding occupation probabilities and 
energies at each point on the energy surfaces shown in Figs. [5] and El the mass parameters 
and the moments of inertia of the collective Hamiltonian are calculated as functions of the 
deformations /3 and 7. The diagonalization of the hamiltonian yields excitation energies and 
collective wave functions that can be used to calculate various observables. The (red) circles 
in Figs. [7] and [8] denote the Q a values and half-lives, respectively, computed for transitions 
0+ s —> 0+ s between eigenstates of the collective Hamiltonian. The differences with respect 
to mean-field values are not large, especially for the heaviest, weakly oblate deformed or 
spherical systems. For the lighter prolate and more deformed nuclei, the differences can be 
as large as the deviations from experimental values. 

Finally, in Figs. [9] and UM we plot the Q a values and half-lives for a chain of odd-even 
and odd-odd superheavy systems, respectively, in comparison with available data 24], |26 |. 



The Q a values are computed as energy differences between mean-field axial RHB minima of 
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FIG. 10. (Color online) Same as in the caption to Fig. [9] but for the a-decay chain of 288 115, and 
for the nucleus 294 11 7. 



parent and daughter nuclei. The minima are determined in the blocking approximation for 
the odd proton (odd-even nuclei), and the odd proton and odd neutron (odd-odd nuclei). The 
equilibrium deformations are determined from the axial energy minima of the corresponding 
even-even systems, and the equilibrium odd-even (odd-odd) configuration is the one that 
minimizes the RHB total energy of the odd-even (odd-odd) system by blocking the odd- 
proton (odd-proton and odd-neutron) Nilsson orbitals. Because one needs to block several 
Nilsson orbitals in order to find the energy minimum, the calculation for odd-even and odd- 
odd superheavy nuclei is restricted to configurations with axial symmetry. Also in this case 
the model reproduces the mass dependence of Q a values and half-lives, and only for the 
heaviest nuclei with Z = 115 and Z = 117 we find significant differences with respect to 
data. The half-lives are calculated with the Viola- Seaborg- type formula of Refs. j^l, I54} ]. 
which takes into account the effect of the odd nucleons by reducing the transition energy 
with respect to Q a by the average excitation energy of the daughter nucleus. This correction 
is necessary when considering ground-state to ground-state Q a values, because the half-life 
is determined by the most probable transition and this occurs between states with the 
same structure (same quantum numbers for odd-even and odd-odd nuclei). We note that 
a somewhat better agreement with data for 293 117 and 294 117 was obtained using the HN 
macro-micro model that includes a more realistic deformation space 52]. From the energy 
surfaces shown in Figs. [5] and EJ one notes that the corresponding even-even systems are soft 
both in (3 and 7 and, therefore, pronounced effects of core polarization can be expected in the 
odd-even and odd-odd nuclei. The detailed structure of these soft nuclei cannot, of course, 
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be reproduced by the simple blocking approximation assuming axially symmetric shapes, as 
used in the present calculation. Nevertheless, the level of agreement with data on the Q Q 
values and half-lives for odd-even and odd-odd nuclei reflects the underlying structure and 
ordering of proton and neutron quasiparticle states predicted by the functional DD-PC1 and 
the separable pairing interaction. 



IV. SUMMARY AND OUTLOOK 

The framework of relativistic nuclear energy density functionals (EDFs) has been applied 
to a study of deformation effects and shapes of superheavy nuclei. The microscopic self- 
consistent calculation is based on the EDF DD-PC1 33j, and a separable pairing interaction, 
used in the relativistic Hartree-Bogoliubov (RHB) model. 

In addition to the equation of state of symmetric and asymmetric nuclear matter, the 
functional DD-PC1 was adjusted exclusively to the experimental masses of 64 axially de- 
formed nuclei in the regions A « 150 — 180 and A w 230 — 250. It is therefore interesting to 
note that, in self-consistent constrained calculations that include axially symmetric, triaxial, 
and reflection-asymmetric shapes, this functional reproduces not only the empirical masses, 
Q a values, and equilibrium quadrupole deformations of actinide nuclei, but also the heights 
of the first and second fission barriers, as well as excitation energies of fission isomers. This 
is an important result as many modern functionals or effective interactions, used in studies 
of superheavy nuclei, are specifically adjusted to data on fission barriers and fission isomers. 

After testing the theoretical framework in the actinide region, we have performed a self- 
consistent RHB calculation of triaxial shapes for two a-decay chains superheavy nuclei, 
starting from 298 120 and 300 120. In both chains the heaviest systems display soft oblate axial 
shapes with minima that extend from the spherical configuration to \^2o\ ~ 0.4 [Z = 120) 
and |/?2o| ~ 0.3 [Z = 118). The intermediate nuclei with Z = 116 are essentially spherical 
but soft both in /3 and 7, whereas prolate deformed mean-field minima develop in the lighter 
systems with Z = 114, Z = 112 and Z = 110. The theoretical Q a values reproduce 
the trend of the data, with the largest difference between theoretical and experimental 
values of less than 1 MeV. Alpha-decay half-lives are calculated using a five-parameter 
phenomenological Viola-Seaborg-type formula, and a good agreement with experiment is 
obtained for both decay chains. The Q a values and half-lives for two chains of odd-even and 
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odd-odd superheavy systems, have been computed using a simple blocking approximation 
in axially symmetric self-consistent RHB calculations. The theoretical values are in rather 
good agreement with the experimental Q a s and half-lives, and only for the heaviest nuclei 
with Z = 115 and Z = 117 one finds significant differences, most probably caused by the 
restriction to axially symmetric shapes. 

For the two chains of even-even superheavy nuclei we have also explicitly considered 
collective correlations related to symmetry restoration and fluctuations in the quadrupole 
collective coordinates. These correlations are implicitly taken into account when adjust- 
ing energy density functionals to binding energies, but their explicit treatment could be 
important in cases when the initial and final states of an a-decay have markedly different 
deformations or shapes. Using a collective quadrupole Hamiltonian based on REDFs, we 
have calculated a-transition energies between ground states of even-even nuclei (0 + — > + 
transitions), rather than between mean-field minima. The resulting Q a values and half-lives 
do not significantly differ from the corresponding mean-field values, except for lighter prolate 
and more deformed nuclei, for which the differences can be as large as the deviations from 
experimental values. 



Together with the recent RMF+BCS study of fission barriers and fission paths of Ref. (23], 
this work has demonstrated the potential of the new class of semi-empirical REDFs for 
studies of shape coexistence and triaxiality in the heaviest nuclear system, including the 
explicit treatment of collective correlations using a microscopic collective Hamiltonian. This 
opens the possibility for a more detailed analysis of this region of SHN, including all presently 
known nuclides with Z = 110 — 118, as well as spectroscopic studies of nuclei with Z > 100. 
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